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I. INTRODUCTION 

Recently, the possibility of the existence of new physics 
above the TeV scale has been considered through the 
introduction of unparticles [HQ. In this scheme one 
admits a hidden sector with a nontrivial infrared fixed 
point A[/, below which scale invariance is explicit. In 
the ultraviolet (UV) regime, at energies above Ajj, the 
hidden sector operator Ojjv of dimension djjv couples 
to the standard model (SM) fields through an opera- 
tor Osm of dimension n via nonrenormalizable interac- 
tions O uv S M/Mfj uv+n ~ A , where M v is the mass of the 
heavy exchanged particle. Below Ajj, the hidden sector 
becomes scale invariant and the operator Ouv mutates 
into an unparticle operator Ojj with nonintcger scaling 
dimension d u . The coupling of field operators can be 
generically written as 
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M d uv +n- 
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The operator Ojj could be a scalar, a vector, a ten- 
sor or even a spinor. Collider signatures I3L |4]as well 
as otherphenomenological aspects [E H, H, [E EE 
El El EE EI EE EE El Elf resulting from this sce- 
nario have been investigated. Also, several cosmologi- 
cal [H HE [H HI and astrophysical [H, [H, [H con- 
straints on unparticle physics have been studied, includ- 
ing bounds arising from stellar equilibrium [26| and black 
hole evolution [27[. The exchange of unparticles gives 
rise to long range forces which deviate from the inverse- 
square law (ISL) for massless particles due to the anoma- 
lous scaling of the unparticle propagator. For example, 
the exchange of scalar (pseudoscalar) unparticles can give 
rise to spin-dependent long range forces, as pointed out 



in Ref. [28?]. Coupling between unparticles and vector or 
axial- vector currents have been investigated in Ref. [23]. 
In Ref. [30( the coupling between unparticles and the 
energy-momentum tensor was studied . Torsion-balance 
experiment results have been used to put limits in these 
interactions [IE [IE [3E HH . In this work we will examine 
the deviations from ISL due to tensor and vector parti- 
cle exchange and the possible constraints that can be 
derived using the bounds on the variation of the gravita- 
tional coupling G, at the time of big bang nucleosynthesis 
(BBN). 



II. UNGRAVITY FROM TENSOR 
UNPARTICLES 

If Ojj is a rank-two tensor it could couple to the stress- 
energy tensor T^ v and its exchange between physical par- 
ticles could lead to a modification of Newtonian gravity. 
Taking the gravitational coupling of the tensor unparticle 
to to be of the form 



(2) 



where M* = Ay (My / A\j) uv , it can be shown that, 
in the nonrelativistic limit and for d u =/= 1, the effective 
gravitational potential of the unparticle exchange has the 
form [H 



V(r) = -G N 



mi 7712 



r 



2d „ 



(3) 



Here Gn should be identified with the Newtonian grav- 
itational constant, Gat = 6.7 x 10~ 39 GeV -2 . The char- 
acteristic length scale Rg for which the ungravity inter- 
actions become significant is defined to be 
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FIG. 1: (color online) Upper bound on Rg as function of d u 
for \ AG/G\ < 0.2, 0.15, 0.1, 0.086, 0.036, 0.01 (from top to 
bottom). 



where Mpi = 1.22 x 10 19 GeV is the Planck mass, and 
C(d u ) is given by 



C(d u ) 



2 T(d u + ±)T(d u - §) 
r 2d --i T(2d u ) 



(5) 



The case d u < 1 leads to forces which fall slower than 
gravity and can be easily ruled out from fifth force ex- 
periments [32| . Hence we will consider only d u > 1 (see, 
however, Ref. [26|). Torsion-balance experiment searches 
for power-law modifications to the ISL have been used to 
constrain the modified potential [3(J. Also, the observed 
perihelion precession of Mercury has been used to test 
these interactions (33j . 

The modification of the gravitational potential can be 
seen as a dependence of the gravitational coupling on r. 
Indeed, the force associated with the potential of Eq. ([3]) 
is 



F(r) = -W(r) = -G(r] 
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with an effective gravitational coupling given by 

Rg 
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In this work we investigate the possible limits on the 
different energy scales (Ajj and AT*) that can be derived 
using the bounds on the variation of the gravitational 
coupling G, 
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at the time of BBN. The effect of a varying G on BBN 
is through the Hubble expansion law, H = a/ a oc \[G 
(a being the scale factor), which determines both the 
neutron-proton density ratio at freeze-out and the effi- 
ciency of 2 H burning abundance when nucleosynthesis 



starts. Given the large statistical and systematic errors 
of the measurements, the typical constraints on the vari- 
ation of G are of the order of a few percent. In particular, 
it is found that 34 1 



AG 

-0.036 < — < 0.086 



(9) 



where the values Y, 



v — 

2-87lo.o2i have been used for 



at 95% confidence level (C.L. 

0. 250 ± 0.003 and 2 H/H 
the 4 He mass fraction and the normalized deuterium 
number density, respectively. The result is dominated 
by the effect on 4 He. Similar results are found by several 
other authors, depending on the experimental values for 
light nuclei adopted in their analysis (see 34] and refer- 
ences therein). In the pessimistic case of a large system- 
atic error on Y p , deuterium can provide a 20% bound [34j . 

In order to infer on the range and strength of the un- 
particle force one needs to know the typical distance r 
between particles interacting during BBN. It turns out 
that during this epoch neutrons and protons have mean 
free paths of the same order of magnitude, A p ~ A„ ~ 
A = 1 m, which in the following we take as the typical 
distance r between nucleons [35l |36||. 

From Eq. © it is now easy to derive upper bounds on 
the characteristic length scale Rg for which the ungrav- 
ity interactions become significant. On Fig. [T] the upper 
bounds on Rq as a function of d u are depicted, consid- 
ering |AG/G| < 0.2, 0.15, 0.1, 0.086, 0.036 and 0.01 . 
On Fig. [5] we show the allowed regions (above the solid 
lines) on the M* — d u (left panel) and M* — Ajj (right 
panel) parameter spaces, which can be derived from 
Eq. (j4|) when combined with the previous constraint. 
On the left panel we present lines for Ajj — 100 GeV, 

1, 10, 100 TeV; and on the right one it has been consid- 
ered d u — 1.1, 1.25, 1.5 and 2. In both cases we have 
fixed AG/G < 0.086. 

These limits should be compared to the ones that 
can be derived from ISL deviations using precision sub- 
millimeter tests. Following [3(3|, we use the data from 
Ref. [13], where the potential 
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is considered, to obtain an upper limit on Rg, by iden- 
tifying d u = (k + l)/2 and R G = /^ /(fe_1) . From that 
upper bound we easily get the constraints on the energy 
scales. Those limits are plotted on Fig. [2] (dashed lines), 
together with the ones obtained from BBN. We interpo- 
lated the limits in Table I of Ref. [13] to obtain bounds 
on Rq as a function of d u . For k < 3 we consider /3k 
as a function of l/(fc — l) 2 , while for larger k we simply 
interpolated the limit linearly. 

We find that these bounds on M* are stronger than the 
ones obtained from BBN. Notice however that, unlike the 
laboratory bounds, our results test ungravity at the early 
universe. 
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FIG. 2: (color online) Allowed region (above the curves) on the M* — d u (left) and M* — A u (right) parameter spaces for 
AG/G < 0.086. On the left it has been taken A v = 100 GeV, 1, 10, 100 TeV [lighter (top) to darker (bottom) curves] and on 
the right d u = 1.1, 1.25, 1.5 and 2 [lighter (top) to darker (bottom) curves] obtained from BBN bounds (red solid lines) and 
ISL violation data (green dashed lines). 



III. LONG RANGE FORCES DUE TO VECTOR 
UNPARTICLES 

Let us now consider the long range forces resulting 
from the coupling of vector unparticles [29]. The po- 
tential for a coupling between a vector unparticle and a 
baryonic (or leptonic) current J M of the form 



A 



A 



j 



(11) 



case we consider the bound AG/G > —0.036. Following 
the line of thought of the previous section we are able 
to derive the limits on A presented in Fig. [3J where the 
solid (dashed) lines represent the upper bounds on this 
coupling as a function of d u and Ajj resulting from BBN 
(torsion-balance experiments) constraints. As in the ten- 
sor exchange case, we find that the bounds arising from 
laboratory searches of putative violations of the ISL are 
more stringent. 



is given by [2£ 
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where are the total number of baryons of the two 
interacting objects and 
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Combined with the gravitational potential we can write 
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with 



Rg = 



1 fXMpi 



A v V u 



C(d u ) 



(14) 



(15) 



where, in order to obtain numerical results, we have made 
the approximation Ai^ ~ mi^/u, u — 931.4 MeV being 
the atomic mass unit. 

The vector unparticle exchange causes a negative vari- 
ation of G, given that the force is repulsive. Hence in this 



IV. CONCLUSIONS 

In this work we have examined the importance of the 
existence of unparticles on BBN yields through the mod- 
ification they introduce in the ISL. We have considered 
tensor and vector unparticle exchange. 

We find that in both cases the BBN bounds are less 
stringent than the laboratory ones searching for viola- 
tions of the ISL. For d u close to unity, the bounds are 
comparable. From the constraints on the variation of G 
during BBN, considering Ajj = 1 TeV and d u = 1.1, 
we find M* > 6.04 x 10 17 GeV and A < 3.54 x 10~ 18 , 
for tensor and vector exchange, respectively, while fifth 
force experiments yield M* > 2.83 x 10 18 GeV and 
A < 1.17 x 10 -18 . We should notice that the range 
Ajj > 1 TeV is the most interesting case to consider since 
it is the one where unparticles may be observed in future 
colliders. 

The difference between BBN and laboratory bounds 
becomes more visible for larger values of d u . For d u = 2, 
we get M* > 15.9 GeV and A < 1.34 x 10" 1 from BBN, 
and M* > 2.36 x 10 5 GeV and A < 1.40 x 10~ 5 from 
torsion-balance experiments. We remark however that 
our bounds concern unparticle effects on the early uni- 
verse and, given the methodological differences, should 
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FIG. 3: (color online) Allowed region (below the curves) on the A — d u (left) and A — A„ (right) parameter spaces for AG/G > 
-0.036. On the left we took Ajy = 100 GeV, 1, 10, 100 TeV [lighter (bottom) to darker (top) curves] and on the right 
d u = 1.1, 1.25, 1.5 and 2 [lighter (bottom) to darker (top) curves] obtained from BBN bounds (red solid lines) and ISL 
violation data (green dashed lines). 
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